Electronic and magnetic properties of the (001) surface of hole-doped manganites 
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The electronic and magnetic properties of ferromagnetic doped manganites are investigated by 
means of model tight-binding and ab initio self-interaction corrected local spin density approxi- 
mation calculations. It is found that the surface alone by breaking the cubic symmetry induces a 
difference in the occupation of the two e g orbitals at the surface. With ab initio calculations we 
found surface localisation of one orbital and hence a change in the Mn valency from four in the 
bulk to three at the sub-surface. Different surface or disordered interface induced localisation of the 
orbitals are considered too with respect to the nature and the strength of the magnetic exchange 
coupling between the surface/interface and the bulk-like region. 



I. INTRODUCTION 

The half-metallic properties of Lai-^Sr^MnOs (x~0.3) 
(LSMO) are of great importance for applications in 
spintronics. The tunnel magnetoresistance (TMR) of 
LSMO/SrTi0 3 (STO)/LSMO junction shows a magne- 
toresistance ratio in excess of 1800% 1 and was, by these 
authors, attributed to half-metallicity. While the mag- 
netic impurities which might diffuse to the insulating 
layer could play an important role in the tunnelling 
process^ through the spin flip effect, antiferromagnons 
at the interface^ due to an antiferromagnetic interlayer 
coupling with the subsurface (or bulk) are found also to 
affect the MR in manganites tunnel junctions. The mag- 
netic properties of these materials are highly sensitive 
to local crystal properties. The extrinsic strain field in- 
duced by lattice mismatch with the substrates or tunnel 
barriers can be sufficient to severely degrade the ferro- 
magnetic order in the surface layers which are critical for 
tunneling^. Although other defects such as segregation 
of a particular species, like Sr in LSMO, at the inter- 
face alters the desired electronic and magnetic proper- 
ties we have not addressed the issue in this work. Good 
TMR is expected if the material is fully polarised and 
half-metallic. Thus the occurrence of an antiferromag- 
netic layer or a localised layer at the surface will be very 
detrimental to tunneling. There is then a strong case 
for demanding to have both half-metallicity and ferro- 
magnetic exchange at the interfaces between manganites 
and insulating barriers. Understanding the spin polari- 
sation at the surface is then of major importance assum- 
ing that growing techniques could fabricate sharp well- 
defined interfaces^ and low diffusion rates of the mag- 
netic ions into the insulating layer. In this paper we 
investigate the conditions required for surface antiferro- 
magnetism and/or localised surface states so as to make 
clear the physical conditions for avoiding them. 

On the basis of local spin density (LSD) band the- 
ory calculations, the origin of half-metallic charac- 
ter of manganese perovskites was discussed in several 
papersiSiLSi2ii2iii These LSD calculations failed to obtain 
a half-metallic state and subsequently the possibility of 
transport half-metallicity was raised by Nadgornyet ali& 



and Mazinjii It could equally well be argued that the 
fascinating electronic and magnetic properties of LSMO, 
including colossal magnetoresistance (CMR), might indi- 
cate that the electronic structure is more complex than 
the standard band theory picture (see reviewsiSiH) and 
might necessitate a better treatment of correlation ef- 
fects. Of particular importance would be to see if these 
correlation effects confirmed the half-metallicity of these 
materials. 

Recently two of us^* described how upon Sr doping of 
LaMn03 (LMO) the Mn valence increases from 3+ to 4+ 
by delocalizing the e s electron. These results therefore 
suggested that, in LSMO, Sr hole doping favours band 
formation instead of localisation. With this Sr doping 
no half- metallic state was obtained in LSMO. Rather, 
the calculations suggested, half-metallicity is the conse- 
quence of remaining local Jahn- Teller distortions from 
the LMO parent material. This did go hand in hand 
with a mixed valence Mn 3+ /Mn 4+ ground state which 
was discovered for Sr concentrations less than 20%^ 

The importance of the local distortions in LSMO does 
suggest that the surface properties of LSMO might be 
different from the bulk properties. This could possibly 
have important consequences for the magnetoelectronic 
transport through an LSMO/STO interface. Surfaces 
of manganese perovskites were studied before: Fillipetti 
and Pickett used a pseudopotential method to study the 
magnetic properties of the surface of CaMn03 (CMO) 15 
and Lai-xCa^MnOa (LCMO)i& in the (001) direction, 
Evarestov et al*£ used the Hartree-Fock approach to 
study the surface (110) in LaMnOa. In particular the 
work of Fillipetti and Picket1ii^i& stressed the impor- 
tance of spin flip processes at the surface for the transport 
properties. 

In this paper we deal with finite slabs of Ri-a, A^MnOa 
(where R and A are trivalent and divalent ions respec- 
tively) to study the surface in a tight-binding model and 
Self-Interaction Correctedi&±£ (SIC) LSDA calculations. 
The first allows us to study larger systems and more com- 
plicated magnetic configuration using a few parameters 
whereas the latter is parameter free and therefore more 
accurate but limited by the number of atoms that can be 
simulated. 
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II. TIGHT-BINDING METHODOLOGY 

The active orbitals in a model calculation on mangan- 
ites are the two degenerate e g orbitals separated by a 
"strong" ligand held from the three low-lying t 2g states. 
As we are concerned with a region of the phase dia- 
gram where most manganites are found to be in the 
FM metallic phase we use the Kondo-lattice type model 
Hamiltonian22i2i using the two e g orbitals: 
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The Hamiltonian consists of the kinetic energy of the 
e g electrons with anisotropic hopping integrals , (7 
and 7' denote the two e g orbitals, i and a index the sites 
and the first neighbours respectively), a Hund coupling 
which favours the alignment of their spins (si) with the 
core- like t 2g moments (Si) and superexchange interaction 
between the classical t 2g spins. 

The transfer integrals between two orbitals e g ( 3z 2_ r 2) 
(orbital 1) and e g ^ x i_ y i^ (orbital 2) on adjacent Mn ions 
are given by: 
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where V£ d is the Slater-Koster parameter 



and e p , ed are the energies of the O 2p and Mn 3d states. 
This parameter takes into account the hybridisation with 
O which does not appear explicitly in the model. 

As shown above an electron in the state e s ( 3z 2_ r 2) can- 
not hop along the z direction whereas its hopping integral 
along x and y is larger than for the e g ^ x 2_ y 2y This fact 
will be important in the determination of the occupancy 
of the two orbitals in the presence of the surface and/or 
an interlayer antiferromagnetic coupling. The strong on- 
site Hund coupling will favour the alignment of neigh- 
bouring core spins via the itinerant e g electrons. Com- 
peting with this tendency is the superexchange which 
acts between core spins on neighbouring sites. This in- 
teraction is responsible for the observed G-AF phase in 
the end members AMn03 where the e g electrons are ab- 
sent and thus only the superexchange operates. The su- 
perexchange also wins over when there are not enough 
carriers to lower the total energy by a gain in kinetic en- 
ergy or in case where hopping is suppressed due to other 
factors as is the case in the presence of a surface as we 
will see below. In order to keep the number of parame- 
ters to a minimum we did not include the Coulomb on- 
site repulsion between e g electrons nor the Jahn- Teller 



coupling 21 . On the other hand we added a shift 21 A of 
the on-site energy for the orbitals at the surface in order 
to take into account the change in energy of the states 
at the surface due chemical shifts and/or strain fields. 
This may be the case at interfaces with grain bound- 
aries or with insulating barriers in tunneling devices as 
explained above. One of the major effects of the surface 
is the occurrence of a charge transfer to or from the bulk 
region inducing a loss of local neutrality and creation of 
electrostatic dipoles. We treat these interactions in the 
Hartree approximation 21 by solving the Poisson equation 
with e = 5. The coupled Schrodingcr-Poisson equations 
are solved self consistently until the relative change in 
energy and charge is less than 0.05%. 



III. RESULTS AND DISCUSSION 
A. Orbital ordering and surface magnetism 

We will consider three possible magnetic configurations 
namely ferromagnetic (FM), A- type antiferromagnetic 
(A-AF) and a configuration where all the moments on 
the inner layers are parallel and the surface one flipped 
(DUUD, where D is for down and U is for up) in or- 
der to look at the interplay between magnetic ordering 
and orbital ordering and the effect of the surface. In 
section III.C we evaluate the energy of a reversed layer. 
We assume an in-plane ferromagnetic ordering so that we 
have one incquivalcnt atom per plane. We show in Fig. ^ 
the occupancies of the e g ^ x 2_ y 2^ and e s ( 3z 2_ r 2) orbitals in 
these three configurations. There is a noticeable correla- 
tion between the type of magnetic and orbital orderings. 
The e g ( x 2_ y 2j is more populated than the e g ( 3z 2_ r 2) on 
those planes which are antiparallel to their neighbours. 
Whereas the two are equally occupied when the coupling 
is FM. The surface layer is an exception however but 
this is not surprising having in mind the anisotropy of 
the transfer integrals defined above. 

The higher occupancy of the e g / x 2_ y 2-\ orbital at the 
surface is explained by the absence of interlayer hop- 
ping for the electrons in this state which do not lose 
kinetic energy in the presence of the surface. Whereas 
the e g ( 3z 2_ r 2) electrons are more sensitive to the pres- 
ence of the surface which limits their hopping and as a 
result this orbital is more occupied in the bulk where the 
levels are broadened than at the surface where the level 
is more localised. As mentioned above the anisotropy of 
the hopping integrals leads to no direct electron transfer 
from e g ( x 2_ y 2} orbitals between planes. Hence in the cur- 
rent model the local density of states (DOS) projected on 
this orbital is independent both of the position in the slab 
and the magnetic orientation of the neighbouring planes. 
On the other hand there is transfer between e 9 ( 3z 2_ r 2) 
orbitals along z. This means that the local e g ( 3z 2_ r 2) 
DOS will be narrowed at the surface because the trans- 
fer is only to one plane instead of two and this is true 
also if there is AF order. We see this effect clearly in 
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FIG. 1: Occupancy of the e s ( 3z 2_ r 2) and e s ( a .2_ H 2) and their 
sum as a function of the distance from the surface for a 21- 
layer TB model of Ro.7Ao.3Mn03. The configurations are in 
order FM, FM with the surface layer flipped (DUUD) and 
A-AF. 



Fig. n Conversely, the decrease in the kinetic energy of 



the i 



g(3z 



r 2-| electrons at the surface results in the weak- 



ening of the double exchange and a tendency to an AF 
coupling between the surface and subsurface layers. The 
higher occupancy of the e g i x 2_ y 2\ orbital at the surface 
makes it more likely that it will want to localise as we 
will see below in the SIC calculations. Inside the bulk- 
like region of the slab the two orbitals are equally likely 
to be occupied as long as the coupling between layers 
is FM. The strong on-site Hund's coupling will always 
act to align the e g electron's spin to the t% g one and in 
order for the system to gain from the kinetic energy of 
the electrons, the core spins ought to be parallel. If this 
is not the case then the hopping is partially suppressed 
and the superexchange wins over. This suppression oc- 
curs between layers and as the only orbital which has a 
finite transfer integral in this direction the e g ^ z 2_ r 2^ will 
be penalised and hence depopulated as can be seen from 
Fig. 2] In Fig. n wnere n0 shift of the surface levels is 



introduced the e 



<?(3z 2 -r 2 ) 



is disfavoured near the surface 



because of the narrowing of the local DOS. The electrons 



would rather go to the wider e g t x 2_ y -2\ DOS. Thus there 
is an orbital order induced at the surface and is present 
in all of the three configurations but is stronger when 
there is a local AFM coupling between the surface and 
subsurface layers as is the case in (b). The effect is even 
bigger in (c) where the AF coupling all through the slab 
causes a narrowing of the e ff ( 32 2_ r 2-) and no significant 
change to the e g r x 2_ y 2y DOS so that we found orbital 
order throughout the film. The total electron density 
shown by the sum, d 3z 2_ r 2 + d x 2_ y 2, is almost unchanged 
at 0.7. This is due to the suppression of charge imbal- 
ance by including the electrostatic interactions solved for 
in the Poisson equation. Without this extra electrostatic 
term the electrons in the absence of a shift of the surface 
levels would transfer to the inner part of the slab raising 
the density above the bulk level independently of the slab 
thickness. We see however from Fig. ^that we get the 
bulk properties for three layers away from the surface. 

We have studied the effect of introducing the shift A 
for both orbitals and for only one of them at the surface. 
As explained earlier, this parameter is used to mimic the 
effect of structural (like lattice mismatch) and chemical 
(like charged surfaces) on the on-site energy level of the 
orbitals. Here we look at the effect on the orbital order- 
ing and in Section III.C we will consider the changes in 
the relative stability of the two solutions FM and DUUD. 
The orbital occupancies are given in Fig|3 as a function 
of the strength A for the two magnetic solutions corre- 
sponding to (a) and (6) of the Fig. ^ We found that 
when both orbitals are shifted by a small A the occu- 
pation of the e g i x 2_ y 2\ remains higher and is explained 
from the simple kinetic energy gains argument explained 
earlier. Increasing A results in a crossover to higher oc- 
cupancy of the e g (3 Z 2_ r 2) orbital but this is explained 
also from the anisotropy of the transfer integrals argu- 
ment. As is depicted schematically in Fig. [3] the LDOS 
of the e g ( x 2_ y 2j orbital is broader than that of e g ( 3z 2_ r 2^ 
so that when the Fermi level is well below the centre of 
the bands the occupation of the first is higher. With 
increasing A the Fermi level lies at the centre and the 
two orbitals are equally filled. Increasing further A will 
favour the occupation of e g ^ z 2 _ r 2j which has higher num- 
ber of states available in a much smaller energy window. 
The crossover occurs for smaller values of A in the FM 
solution than in the DUUD case, in agreement with our 
earlier finding that the occupation of e g ^ z 2_ r 2^ favours 
FM coupling. Shifting one orbital only will favour its 
occupation in both solutions and for values of A larger 
than 2to the other orbital is completely depleted. 



B. Ab initio study of charge and orbital ordering 

As mentioned above we report also on results using 
the SIC-LSDA to study the surface of a a representative 
system, Lao.7Sro.3Mn03. This method has already been 
used with success in studying the bulk properties of this 
material-- 4 . It allows for a parameter-free total energy 
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2: Evolution of the occupancies of the e 



g 0c 3 -y 3 ) orbitals at the surface versus the shift A. 
figurations are considered: FM(a) and DUUD (b). 
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FIG. 4: LDOS, in the rigid band model, for a 5 units su- 
percell (Lao.7Sro.3Mn03)5. The Mn in layers 2 and 10 have 
localised 3i2 9 + le g ^ x 2_ y 2- j electrons and the other Mn have 
localised 3t2 9 electrons only. The left-hand-side picture shows 
the supercell. 




FIG. 3: Schematic representation of the evolution of the 
occupancy of the e ff ( 3z 2_,.2) (narrow) and e g ( x 2_ y 2) (broad) 
orbitals at the surface when shifted by an amount A with 
respect to their common bulk level. The area below the two 
curves is the same. The horizontal line represents the Fermi 
energy. When the Fermi level is above the centre of the two 
bands the occupation of e 9 ( 3z 2_,.2) becomes higher 



minimisation with respect to the localised/itinerant state 
of the Mn d electrons in our case. The number of elec- 
trons allowed for band formation is found by comparing 
total energies in the two configurations where the elec- 
tron is itinerant and where it is localised. The valency 
of the Mn ions is then found by subtracting the localised 
electrons from the total number of valence electrons. It 
has been applied successfully to systems where there is 
strong tendency to localisation of the valence electrons 



that could not be accounted for using the conventional 
LSDA functionated 

Here the focus is put on the changes brought about by 
the surface on the charge and orbital orderings. The va- 
lence of Mn in LSMO was calculated to be tetravalenfei^ 
in the bulk material. However we find here that this 
valence is reduced to trivalent when the Mn is on the 
subsurface layer. We report on calculations in which we 
include a virtual La/Sr atom to account for the mixed va- 
lence. All the layers are either MnC>2 or Lao.7Sro.3O. This 
is a type of rigid band model. First we studied a super- 
cell consisting of four layers of MnC>2 with three and four 
layers of Lao.7Sro.3O. The second system is considered 
in order to check the effect of the stoichiometry which 
is not respected in the system with three Lao.7Sro.3O 
layers only. This non-stoichiometric system is symmet- 
ric and the calculations are much less involved than in 
the stoichiometric but non-symmetric case. The slabs 
are separated by seven and six layers of empty spheres in 
the symmetric and non-symmetric case respectively. The 
non-symmetric system has two surface terminations, i.e 
Mn02 and Lao.7Sro.3O whereas in the symmetric case 
the termination is Mn02- 

In the LSDA calculations we found an energy differ- 
ence of 8.04 mRy/(Mn02 layer) between the ground- 
state FM configuration and the configuration where the 
surface moment is antiparallel to the bulk ones (DUUD) 
in the symmetric case. For the non-symmetric case the 
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ground state is also FM and the difference in energy 
with the flipped surface moment configuration is of 5.59 
mRy/(Mn02 layer). The surface and subsurface Mn 
magnetic moments in the ground state are of 3.26 and 
3.16/is in the symmetric system and of 3.28 and 3.12/1^ 
in the non-symmetric system. The results are indeed in 
good agreement. We then studied different orbital local- 
isation scenarios. These are the localisation of the 3^2g 
orbitals all through the slab and localisation of an extra 
e g orbital at the surface, which gives three possible sce- 
narios for each of the magnetic configurations. In both 
systems the ground state is the FM phase with the Zti g 
orbitals only localised on all the Mn ions. The FM con- 
figuration with an extra e g ( x 2_ y i^ localised at the surface 
is 6.67 mRy/(Mn02 layer) higher in the symmetric case 
and is of 6.21 mRy/(Mn02 layer) in the non-symmetric 
case. The surface and subsurface Mn magnetic moments 
are of 3.30 and 3.31/Lts in the symmetric system and of 
3.32 and 3.27 fis in the stoichiometric system. 

Having confirmed that the stoichiometry has negligi- 
ble effect on the overall relative stability of different or- 
bital and magnetic configurations we applied this method 
to symmetric five-Mn02 layers supercell of LSMO sur- 
face which included four unit cells of empty spheres (see 
Fig. 0J . In this case the MnC>2 layer is at the sub-surface 
rather than the surface which we studied previously. The 
TB model does not include the Lao.7Sro.3O layers and 
therefore cannot differentiate between these cases. How- 
ever a phenomenological shift A of the on-site energies 
at the surface could capture this. We found that termi- 
nations by the Lao.7Sro.3O layer leads to the localisation 
of one more electron on the Mn atoms in the Mn02 layer 
under surface. The ground state configuration has lo- 
calised 3<2g + l e g(x 2 -j/ 2 ) electrons under the surface and 
3^2g electrons on the manganeses in the bulk in agree- 
ment with the model calculations when the shift A is ap- 
plied to the surface e g r x 2_ y 2^ orbital (see Fig. EJ. Since 
the Mn02 layer is not the termination, the e g ^ z 2_ r 2^ 
could not be said to be favoured as was the case in ear- 
lier calculations^ 6 *^. There is however an electrostatic 
interaction with the surface 2 ~ ion which means that 
the energy level of this orbital should increase in compar- 
ison to the case when the Mn02 layer is at the surface. 
Whereas this is not true for the e g r x 2_ y 2\ orbital which 
is less sensitive to the presence of the extra Lao.7Sro.3O 
layer. The situation can then be modelled by shifting the 
on-site e g i x 2_ y 2-\ level downwards which is equivalent to 
shifting the other level upwards. The calculated ground 
state configuration has the total energy of 37 mRy lower 
than the system with all Mn 4+ configured manganese. 
From the LDOS of the first layer we can see small contri- 
butions of electronic states from the vacuum region. The 
Lao.7Sro.3O surface layer becomes insulating with a band 
gap of about 1 eV. The LDOS for the first Mn0 2 layer 
with Mn 3+ valence shows a nearly half-metallic charac- 
ter with a nearby pseudo-gap of 1.7 eV. This shows that 
localisation of the e g ^ x 2_ y 2^ electron, forced by the sur- 
face, leads to near half-metallic properties at the LSMO 




FIG. 5: Energy differences between the two configurations 
where the surface-subsurface coupling is AFM or FM in units 
of the hopping parameter to per surface Mn ion versus the 
amount A of the shift added to one or both surface levels. 
The Hund's coupling parameter is Jh = 6to- 



surface as was mentioned beforei*. The change of the 
symmetry of the localised electron to e ff ( 3z 2_ r 2) increases 
the total energy by about 23 mRy. This is a substan- 
tial energy, 62% of the energy needed for delocalizing 
this e g electron. Additionally, the rotation of e g ( 3z 2_ r 2) 
orbitals by 90° into the x — y plane makes this configu- 
ration unfavourable by only 10 mRy in comparison with 
the ground state configuration of a localised e g ( x 2_ y 2) 
electron. Localisation of e g electrons on the other Mn 
atoms, inside the bulk, increases the total energy of the 
system. The LDOS for the next LaO layer shows metallic 
character with a small number of electrons at the Fermi 
level but each of the Mn02 layers including the Mn 4+ has 
a clear metallic character. The LDOS for the layers below 
these are similar to the ones for the bulk. We noted that 
decreasing the volume of elementary cell increases the dif- 
ferences between the energies described above. However, 
increasing the lattice parameter by 2% decreases the dif- 
ferences between energies states ten times. The magnetic 
order at the surface changes through spin flip processes 
from the bulk to a local anti- ferromagnetic arrangement. 



C. Energetics 

Depending on the surface/interface termination several 
scenarios are possible which will allow for a particular 
orbital to be favoured. It has been argued that a ter- 
mination with Mn02 plane will favour the occupation of 
the e s (3 Z 2_ r 2) orbital-i which has its lobe oriented toward 
the missing oxygen ion and hence has lower electrostatic 
energy than the e g ( x 2_ y 2-^ which still sees the oxygens 
present in the plane. A phase diagram in the parame- 
ter space has been obtained^ where it is shown clearly 
that a large shift A of the e g ( 3z 2_ r 2) with respect to 
e g(x 2 -y 2 ) an( i the bulk levels will favour an inplane anti- 
ferromagnetism at the surface and a canted configuration 
with respect to the "bulk" . First principles calculations 16 
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on the other hand found that the surface-subsurface ex- 
change coupling in La^Cax^MnC^ is FM independently 
of the "bulk" magnetic configuration. The authors^ ar- 
gued that this is due to the dangling e 9 ( 3z 2_ r 2) bond 
again which favours this ferromagnetism. Whereas in the 
present ab initio study we have found that it is rather the 
e g(x 2 -y 2 ) that is localised as the Mn02 layer is the sub- 
surface rather than the terminating layer. We present 
below results concerning the competition between FM 
and AFM surface-subsurface coupling as a function of the 
symmetry of the localised orbital and the amount A by 
which the levels are shifted. We considered shifts of both 
levels simultaneously and of one orbital at a time and 
calculated the band energies (kinetic plus Hartree) when 
the surface-subsurface coupling is either FM or AFM. 
The superexchange energy difference between them is of 
2Jaf per surface Mn ion. The results are shown on Fig. [5] 
where the energies are given in units of the hopping in- 
tegral to and the Hund's coupling constant is taken 
equal to &t . As can be seen from the figure while local- 
ising both orbitals does not change the relative energies, 
localising the e g ( x 2 -y 2 ) wm affect strongly the surface- 
subsurface ferromagnetism. This is due to the depletion 
of the other orbital on the surface and hence the weak- 
ening of the double exchange mechanism mediated by 
e g(3z 2 -r 2 ) electrons hopping between the layers. In the 
case where both orbitals are shifted it is the reduction in 
the subsurface occupation of the e g ( 3z 2_ r 2) orbital which 
limits the gain in kinetic energy that would result from 
the higher occupation of the surface e g ( 3z 2_ r 2) orbital. 
The relative stability of the two configurations remains 
unaltered as a result. Shifting the e g ( 3z 2_ r 2) will enhance 
the FM coupling with the bulk as found in the ab ini- 
tio study of the MnCVterminated system. But this is 
true only if the e g ^ x 2_ y 2j is not strongly disfavoured as 
found in the previous model calculations^. If the occu- 
pation of the latter orbital is low at the surface in-plane 
anti-parallel coupling of the spins would result as a con- 
sequence of the weakening of the double exchange at the 
surface which is mediated mostly by e g ( x 2_ y 2j electrons. 

IV. CONCLUSION 

We have considered what light our calculations have 
shed on the ideal hole-doped-manganite-insulator inter- 
face such that tunneling magnetoresistance (TMR) is op- 
timal. We have studied changes that are induced by the 
lack of cubic symmetry at the surface as well as different 



chemical environments which favour the formation of lo- 
calised states, through a realistic double exchange model 
and first principles calculations. In the model calcula- 
tion we have taken account of the different scenarios by 
adding a shift A to the surface on-site energy of the or- 
bitals. If a surface/ tunnel barrier has net positive charge 
A will be negative for both orbitals. If this is too large 
we get localisation which is bad for tunneling. Equally a 
strong negative charged termination is also bad because a 
positive A will deplete both orbitals leading to magnetic 
disorder at the surface. If the in-plane lattice constant of 
the barrier is smaller than that of the manganite crystal 
there will be a strain field which gives a negative A 3z 2_ r 2 
tending to favour the e g ( 3z 2_ r 2) orbital. For small values 
of A 3z 2_ r ,2 surface ferromagnetism is enhanced. However 
large values of this strain will deplete the e g ( x 2_ g 2) orbital 
and favour in-plane antiferromagnetism at the surface 
MnC>2 layer which is detrimental to TMR. On the other 
hand strains favouring e g ( x 2 -y 2 ) always suppresses ferro- 
magnetic ordering and hence TMR. In the absence of any 
chemical shift of the atomic levels at the surface only a 
small amount of charge is transfered to the bulk because 
of the resulting electrostatic forces tendency to restore lo- 
cal charge neutrality. Shifting both levels or one of them 
however will result in charge transfer to the surface for 
large values of A. This leads to the formation of Mn 3+ 
at the surface. These findings are confirmed by the more 
accurate ab initio SIC-LDA calculation on a model sys- 
tem Lao.7Sro. 3 Mn0 3 . In these calculations we found no 
localisation of the surface orbitals when the slab is termi- 
nated by a Mn02 layer and the coupling is ferromagnetic. 
We studied the LaSrO-terminated system were the sub- 
surface e s ( 3z 2_ r 2) still sees an oxygen ion on the surface 
and interacts strongly with it. This orbital is then dis- 
favoured and we found that the e g t x 2_ y 2\ is localised at 
the subsurface layer changing the valency of the Mn ion 
from tetravalent in the bulk to trivalent at the surface. 
The magnetic coupling then becomes antiferromagnetic 
as would be expected from the correlation between or- 
bital and magnetic ordering. However if the interlayer 
distance is increased at the surface, which amounts to 
favouring the e g ( 3z 2_ r 2), we found that the energy differ- 
ences between competing orderings decrease significantly. 
This relates the ab initio results to the model ones. In 
summary we predict that the best TMR will come from 
tunnel barriers that are neutral or weakly positive. There 
should be a minimal surface strain although a small in- 
plane compressive strain favouring the e g ( 3z 2_ r 2) is actu- 
ally beneficial. 
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